Enhancing synchronization by direct ionahty in complex networks 
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We proposed a method called residual edge-betweenness gradient (REBG) to enhance synchroniz- 
abihty of networks by assignment of link direction while keeping network topology and link weight 
unchanged. Direction assignment has been shown to improve the synchronizability of undirected 
networks in general, but we find that in some cases incommunicable components emerge and net- 
works fail to synchronize. We show that the REBG method can effectively avoid the synchronization 
failure {R = ^2/^^ = 0) which occurs in the residual degree gradient (RDG) method proposed in 
Phys. Rev. Lett. 103, 228702 (2009). Further experiments show that REBG method enhance 
synchronizability in networks with community structure as compared with the RDG method. 
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Synchronization is an important phenomenon in var- 
ious fields including biology, physics, engineering, and 
even sociology 1]. In particular, synchronization in com- 
plex networks has been intensively studied in the past 
decade P-[Toj . One important objective in these studies is 
to enhance the synchronizability [TwlOi], i.e., the ability to 
coordinate oscillators in synchronization. Many related 
methods based on node properties and link weight have 
been proposed. For example, some researchers took into 
account the degree centrality 0, [1] , the betweenness Q 
and the node age [l^ [ll], and assigned weight on links 
to enhance synchronizability. Nishikawa and Motter pro- 
posed to assign zero weight to particular links which lead 
to an oriented tree with normalized input strength and 
no directed loop [H, [l^. Moreover, they proved that 
tree is the optimal structure on which synchronizability 
is highest. Recently, the shortest oriented spanning tree 
is shown to be the optimal structure for both synchroniz- 
ability and the convergence time [13, [ill- However, these 
methods are all based on assigning weights to links, where 
the influences of link direction on synchronization have 
not been intensively considered [l[. 

How to improve synchronization in directed network 
is still an unsolved problem. Though previous works 
suggest that hierarchical structure and the absence of 
feed-back loop can enhance synchronizability, the under- 
lying mechanism is not clear. On the other hand, the 
directionality plays a significant role in the dynamic of 
networks [16l - [l9| . With the understanding of relations 
between link direction and synchronization, a lot of ap- 
plications can be made [l|. For example, simply regulat- 
ing the direction of the phase signal of alternating cur- 
rent can facilitate phase match in power grids without 
additional construction cost in the topology. With this 
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in mind, the authors in Ref. [20|] proposed the residual 
degree gradient (RDG) method to enhance directed net- 
works' synchronizability by assigning only the direction 
of links without changing the entire topology and link 
weights. They also claimed that RDG method can en- 
hance the network synchronization, contrary to the ran- 
domly assigned directional method. 

However, instead of enhancing synchronizability, we 
find that the RDG method results in incommunica- 
ble components in some particular cases, which lead to 
graphs incapable of complete synchronization. Incom- 
municable components of a network correspond to the 
components on which information cannot be transmitted 
from one to the other in either direction. Under these cir- 
cumstance, the networks can never reach complete syn- 
chronized state. In this paper, we propose the so-called 
residual edge-betweenness gradient (REBG) method to 
resolve the problem of incommunicable components in 
the RDG method. By evaluating the betweenness on 
all edges, we devise an algorithm to assign link direc- 
tion. The effectiveness of the algorithm lies in the use 
of edge betweenness, which embeds global information, 
as compared to the node degree which refiect only local 
information. We find that REBG is effective in enhanc- 
ing synchronizability without leading to incommunicable 
components in most network topologies. 

To begin our analysis, we make use of the Master 
Stability Function which allows us to use the eigenra- 
tio R = Amin/Amax of the Laplaciau matrix to represent 
the synchronizability of a network, with Amin and Amax 
denote respectively the smallest nonzero eigenvalue and 
the largest eigenvalue f2ll - [23| . Specifically, < i? < 1, 
the larger the value of R, the stronger the synchronizabil- 
ity of networks. According to Ref. [l3, HO]) we use the 
real part of eigenvalues from the Laplaciau matrix to in- 
vestigate the propensity for synchronization in directed 
networks. However, instead of setting R = ^^in/^max 
(superscript r denotes the real part of complex number). 



we set R = Aj/A^, with Aj to be the second smahest 
eigenvalue. It is because Aj = in cases when incom- 
municable components emerge, which is not reflected by 
AJ'jjjjj. In the subsequent analysis, we use R to character- 
ize synchronizability. 

As the known fact, when there is isolated node or com- 
munity in an undirected network, the network can never 
reach complete synchronization since there is no informa- 
tion flow between the isolated components. Under this 
circumstances, the synchronizability index R = X2/XN — 
0. We call this phenomenon (R — 0) synchronization fail- 
ure. Generally speaking, synchronization failure is more 
likely in directed networks as isolated components are 
not necessary. We first denote the cut-vertex as the only 
node through which two or more components commu- 
nicate with each other. In directed networks, the syn- 
chronization failure happens whenever cut-vertex have 
only incoming links. In this case, there is absolutely no 
communication between those components and the cut- 
vertex becomes an information sink. Hence, synchroniza- 
tion cannot be achieved among the incommunicable com- 
ponents. 

In order to examine synchronization failure, we first de- 
scribe briefly the RDG method [20l|. The RDG method is 
proposed to enhance the synchronizability of undirected 
networks by simply assigning the link direction. They 
refer links without yet assignment of direction as resid- 
ual edges, and the number of connected residual edges 
as the residual degree of a node. In each step, they se- 
lect the node with the smallest residual degree and set 
a maximum of \{k)/2] of its residual edges pointing to 
it[22|, where (k) is the average degree in the original net- 
work. Once a node has been selected, it will not be cho- 
sen again. Nodes have not yet been selected are called 
residual nodes. In addition, there is a directionality a to 
control the final fraction of links with direction assign- 
ment. When a = 0, all links remain undirected. When 
a = 1, all links are assigned direction. The RDG assign- 
ments are finished when there is no residual node left in 
the network. 

The RDG method may lead to directed graph with syn- 
chronization failure, i.e., cut- vertex with only incoming 
links. A simple example is given in Fig. 1(a). According 
to the rule of RDG, node 1 (fc = 2) will be selected first 
and the two remaining communities will be left incom- 
municable. As we have discussed, this RDG network can 
never reach complete synchronized state. 

According to [13], RDG was supposed to provide a 
practical way for the synchronization improvement in 
power grid networks, which have exponential degree dis- 
tribution [25]. So we first study the RDG method in 
random exponential networks [26| . Random exponen- 
tial networks in our context refer to random networks 
with exponential degree distribution. We also tested the 
RDG method in random scale-free networks since power- 
law degree distributions are widely observed in empiri- 
cal data 26]. The degree distribution are respectively 
given by P{k) ~ e"^*^ and P{k) ~ For each /3 




FIG. 1: (a) A simple example of RDG network with resultant 
i? = and original R = 0.039 (b) The REBG network from 
the same original network with 9 = 1, and resultant R = 0.4. 
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FIG. 2: (Color online) The synchronization failure rate as a 
function of a by the RDG method in (a) random exponential 
networks (P(fc) ~ e""*, = 500, fc^in = 2, /3 = 1) and (b) 
random scale-free networks {P{k) ~ A'^ = 500, femin = 2, 
7 = 5). Given a = 1, the failure rate as a function of /3 in (c) 
random exponential network and the failure rate as a function 
of 7 in (d) random scale-free networks when RDG and REBG 
(61 = 0.01, = 0.2 and 6 = 1) methods are used. 



and 7, we tested the RDG method on 100 network re- 
alizations. To examine synchronization failure, we make 
sure that there is no isolated nodes or communities in the 
original networks which masks the effect of incommuni- 
cable components. We find that many resultant RDG 
networks are with R = 0. The failure rate, i.e. the frac- 
tion of realization incapable of complete synchronization, 
is reported in Fig. 2. As shown in Fig. 2(a) and (b), 
when a increase (i.e., more links are assigned direction), 
failure rate increases for both exponential and power- 
law networks. These results imply that synchronization 
failure is common in both power-grid-like and scale-free 
networks, given RDG is used as the method for direc- 
tion assignment. We further tested the case of a = 1, 
which leads to the largest synchronization improvement 
as reported in [20]. By varying /3 and 7, we find that 
failure is high when /3 and 7 are large. We also remark 
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FIG. 3: (Color online) The difference D — -Rrebg — -Rrdg 
in (a) random exponential networks and (b) random scale- 
free networks with A'^ = 500 and fcmin = 2. The results are 
obtained by averaging 10 independent realizations. 
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FIG. 4: (Color online) The fraction of realizations with di- 
rected loop as obtained by the REBG and RDG method in (a) 
random exponential networks and (b) random scale-free net- 
works. The results are obtained by average 100 realizations. 
The abrupt jumps in random exponential networks come from 
the discontinuity in the assignment constraint [(fc)/2]. 

N . All these suggest that synchronization failure is not 
a rare phenomenon and therefore a new method of direc- 
tion assignment is required to prevent failure. 

We thus introduce the residual edge-betweenness gra- 
dient (REBG) method to solve the problem of synchro- 
nization failure. Instead of the node degree, we take the 
edge-betweenness into account. First of all, we define Si 
for node i as 

N 

Si — ^ ^ ^ijhj ; (1) 
J = l 

where = 1 when there exists an undirected link be- 
tween i and j and otherwise 0. Uj is the betweenness of 
the link between i and j evaluated on the original undi- 
rected networks, subjected to a power with < 9 < 1. 
To assign link direction, we select the node with the 
smallest residual Si in each step and assign an incom- 
ing direction for a maximum of \{k)/2\ of its residual 
links. As more directed links are assigned, the residual 
Si has to be updated at every step. If there are multiple 
nodes of the smallest s^, we choose the node with the 
smallest initial Si first. The REBG method stops when 
there is no residual node left in the network. We remark 
that when the parameter 9 — 0, Si — ki and the REBG 
reduces to the RDG method. When 9 > 0, the REBG 
method contains the global information delivered from 
the edge betweenness. 

An shown in Fig. 1 (b), the REBG method can effec- 



tively avoid the problem of synchronization failure. In 
this simple network, whenever node 1 is chosen before 
node 3 and 7, synchronization failure occurs. We note 
that node 1 is the cut-vertex connecting the two com- 
munities, its edges are of high betweenness. To see how 
failure is avoided, we evaluate the initial Si before any 
direction assignment, as given by si = 20^ -I- 20^, Sx = 
l0^ie^Q0 foj. a; = 2, 4, 5, 6, 8, 9 and Sy = 6^ -h 6^ -f 6^-^20^ 
for y = 3,7. By tracing all the possibilities of subsequent 
direction assignment, we find that synchronization fail- 
ure does not occur provided that node 1 is not selected at 
the first assignment. In this case, si > Sx which implies 
9 > 0.18. We denote this value as 9c which marks the 
value of 9 from which failure ceases. We remark that the 
value of 9c is different for different topology. 

We then examine the failure rate in the random ex- 
ponential networks and the random scale-free networks. 
As shown in Fig. 2(c) and (d), failure rate vanishes for 
both networks when 9 = 0.01, 0.2 and 1. These results 
suggest 9c « 0, i.e., failure ceases when edge betweenness 
is considered in direction assignment. It implies that 9 
in Eq. (1) should be positive to make complete synchro- 
nized state possible. 

Finally, we compare the synchronizability index R be- 
tween REBG and RDG methods. Both random expo- 
nential networks and random scale-free networks are ex- 
amined. For better illustration, we report the difference 
D — i?REBG ~ .Rrdg as a function of 9 and /3 in Fig. 
3(a), and 9 and 7 in Fig. 3(b). The positive D shows 
that the REBG method results in higher synchronizabil- 
ity as compared to the RDG method. This enhancement 
is mainly resulted from the prevention of the synchroniza- 
tion failure. One may question the validity of indicating 
synchronizability by R with only real part of eigenvalues. 
We show in Fig. 4 that for the regime when /? and 7 is 
large, the resultant networks are free of directed loops, 
which support the validity of positive D (as obtained by 
R) in this regime. These results also explain the negative 
D observed with intermediate value of /3 and 7. From the 
lines oi 9 = and 6* = 1, we see that the number of net- 
works with directed loops is higher when 9=1, which 
hinders synchronization and lead to unfavorable results 
from REBG. We further show that the REBG method 
with 6* 0.2 are similar to the RDG in terms of the 
fraction of loopy realizations, suggesting 9 = 0.2 is an 
effective value for direction assignment, and at the same 
time avoids synchronization failure. Hence, one can use 
the REBG method with small 9 to enhance synchroniz- 
ability effectively. 

We further compare the RDG and the REBG methods 
in graphs where failure does not occur. Here, we study 
the two methods in Girvan Newman benchmark (GN- 
benchmark) network which consists of 128 nodes and 
is divided into 4 communities [l^j. In GN-benchmark, 
Winter + ^outor = 16, where /cinter is the average node de- 
gree in each community and fcouter is the average node 
degree among different communities. We show in Fig. 5 
the results of /coutcr = 1 from which the GN-benchmark 
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FIG. 5: (Color online) The synchronizability of the REBG 
method and the RDG method in GN-benchmark(fcoutcr = 1) 
as measured by (a) A2, (b) and (c) R — X^/^^n ^ func- 
tion of 9. (d) The order parameter r{t) of Kuramoto model 
(cr = 10) on the RDG networks and REBG networks(6l = 1). 
The results are obtained by averaging 100 independent real- 
izations. 



networks are highly clustered. We can see from Fig. 
5(a) that the synchronization enhancement comes mainly 
from A2. Moreover, we tested the Kuramoto model 
on the resultant RDG networks and REBG networks. 
The oscillator on node i of the networks is described 
by 9i = (jJi + aJ2i Aij sin(9j — 9i) in which A is the 



adjacency matrix and Aij denote the information flow 
from j to i, and the collective phase synchronization 
can be investigated by the order parameter defined as 
r{t) = (l(Ej^ie*^'^'V^)l)- From Fig. 5(d), it is ob- 
vious that r{t) of the REBG networks converges faster 
than that of the RDG networks. These results show that 
the REBG method lead to greater improvement in syn- 
chronizability than the RDG method in highly clustered 
networks, despite the absence of failure in RDG networks. 

In summary, we introduced the residual edge between- 
ness gradient (REBG) method for direction assignment, 
which overcomes the problem of emergence of incom- 
municable components in the residual degree gradient 
(RDG) method. The effectiveness of our method lies in 
the use of edge betweenness, which reflects global net- 
work information when compared to the node degree in 
RDG. Further tests of REBG and RDG in highly clus- 
tered networks show that REBG can lead to greater syn- 
chronizability improvement in networks despite the ab- 
sence of failure problem. For daily applications, such 
incommunicable components brings huge loss to the elec- 
trical systems when power grids are unable to reach com- 
plete synchronization. Hence, the REBG method is ef- 
fective in improving synchronizability which may lead to 
wide applications. 
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